HOW MODEL SETS CAN BE DETERMINED BY THEIR 
TWO-POINT AND THREE-POINT CORRELATIONS 



XINGHUA DENG AND ROBERT V. MOODY 

Abstract. We show that real model sets with real internal spaces are 
determined, up to translation and changes of density 0, by their 2- and 
3-point correlations. We also show that there exist pairs of real (even 
ID) aperiodic model sets with internal spaces that are products of real 
spaces and finite cyclic groups whose 2- and 3-point correlations are 
identical but which are not related by either translation or inversion of 
their windows. All these examples are pure point diffractive. 

Placed in the context of ergodic uniformly discrete point processes, 
the result is that real point processes of model sets based on real internal 
windows are determined by their second and third moments. 



1. Introduction 

An enduring problem of crystallography is the inference of the internal 
structure of a crystal from physically measurable quantities, notably diffrac- 
tion. Perfect knowledge of the diffraction is equivalent to perfect knowledge 
of the 2-point correlation of the structure. However, even perfect knowl- 
edge of the diffraction, or 2-point correlation, is insufficient to pin down the 
structure of a crystal, with counterexamples going back to L. Patterson [22] 
(see for a good source of information on this subject). 

Quasicrystals present the same problem, but are even more difficult. 
Based on the theory of the covariogram, Baake and Grimm [H [11] have 
given examples of model sets (or cut and project sets, as they are often 
called), which are intrinsically different but have the same diffraction. 

One can ask whether knowledge of additional higher point correlations, 
notably the 3-point correlation measure, could provide sufficient information 
determine the structure. But again, even for periodic structures, there are 
counterexamples [12]. The main result of this paper (Thm. [3]) however, is a 
positive result: in the context of regular model sets with real internal spaces, 
the 2- and 3-point correlations do determine the model set (up to translation 
and to modification by sets of density 0). 

The possibility that in pure point diffractive sets (of which regular model 
are good examples) the 2- and 3-point correlations would suffice to know 
all the higher correlations was suggested by D. Mermin in a very interesting 
paper [18j on a new approach to handling symmetry for quasicrystals. His 
ideas are based on the Landau approach to second order phase transitions 
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and involve rather informal manipulation of quantities which, as the author 
recognizes, cannot be mathematically justified as given. To quote from 
that paper: ... this informal Fourier space argument that the identity 
of all second and third order correlations implies the identity of all higher 
correlations is disarmingly trivial. I would very much like to learn of a 
comparably simple informal argument or an instructive counterexample in 
position space. 

In spite of Thm. [31 Mermin's suggestion is not correct in general. This 
can already be seen in the periodic case from the results of [12], which we 
illustrate here in ^7.11 in the form of periodic model sets. In ^7.21 we offer 
an example of a pair of aperiodic model sets, based on an internal space 
which is the direct product of a real line and a finite cyclic group, for which 
the 2- and 3-point correlations are identical but the point sets themselves 
cannot be transformed into each other either by translation or inversion of 
their windows. 

The situation with more general internal spaces is, no doubt, difficult. 
Even the simple case of a product of a real space and a finite group alluded 
to above is quite complicated. We touch on it here in ^ Although, as 
we have just pointed out, counterexamples can occur in this case, in other 
instances we can obtain a positive result. For instance, it is easy to see 
that the vertices of rhombic Penrose tilings, which are model sets based on 
X Z/5Z, are determined by their 2-and 3-point correlations. 

Aperiodic sets are often studied in the context of dynamical systems 
and/or stochastic point processes. This approach was pioneered in [23] 
and has been used extensively both in mathematical and physical models, 
[HiiniE]- This applies particularly to regular model sets (see ^for defini- 
tions) where things can be stated much more precisely, [131 1251 [2]. Instead 
of a single model set one considers its hull, namely the set of all uniformly 
discrete points sets that are in the closure of its translation orbit (see ^for 
more details). This hull is then uniquely ergodic, and we may view it as 
describing a uniformly discrete ergodic point process. 

Any such uniformly discrete point process (with a common lower bound 
on the distance between closest points) is characterized by knowledge of 
its entire set of moments (second, third, etc.) [7]. Knowledge of the kih. 
moment is equivalent to knowledge of the /c-point correlation. Thus our 
result about model sets with real internal spaces can be rewritten (Thm. H]) 
as the statement that for them only finitely many (namely the second and 
third) of these infinitely many moments are needed. 

The main idea behind proving Thm. [3] is to transfer the correlation prob- 
lem to the internal space of the cut and project scheme defining the model 
sets. There the correlations are directly related to what have been called 
in [13] the A:-deck functions, or the covariogram in the case k = 2. These 
become tractable after Fourier transformation, a fact that has been discov- 
ered several times before, see for example |12l I14j . The main obstacle is 
dealing with the set E of zeros (the extinctions in the diffraction) of the 
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Fourier transform of the characteristic function of the window of the model 
set. Here we offer Prop. [2] which we have not seen exphcitly stated in the 
hterature, which allows us to proceed as long as E has no interior points. 

Finally in ^we point out jlT] which offers a different approach via spec- 
tral theory to determining uniformly discrete ergodic pure point diffractive 
point processes by means of their moments. However, ultimately it too re- 
turns to similar problems about extinctions. For a short survey covering 
this and material on point processes, see |21] . 

2. Model sets 

We work in M*^. The usual Lebesgue measure will be denoted by £. The 
open cube of side length R centred at is denoted by Cr, so l{Cji) = R^. 
A cut and project scheme for is a triple S = {W^,H,C) consisting of a 
compactly generated locally compact Abelian group H and a lattice C C 
M'^ X H for which the projection mappings tti and tt2 from x H onto M*^ 
and H are injective and have dense image respectively: 

(1) U 
L ^ C 

Then L := iri{C) is isomorphic as a group to C (though it is rarely a 
discrete subgroup of ffi*^) and we have the composite mapping (•)* : L — > H 
with dense image defined by 7^2 ° {'^i\c)~^- 

The statement that £ is a lattice is equivalent to saying that it is a 
discrete subgroup of M°' x H and that the quotient group T := (M'^ x H)/C 
is compact. Notice that ii H = {0} then L is a lattice in and we are 
back in the situation of normal periodic crystallography. Thus the theory 
of model sets is a generalization of the theory of periodic sets and includes 
them as special cases. 

We let 9h be a Haar measure on H, scaled so it gives measure 1 to any 
fundamental region of £ in the space M'' x H under the product measure 
If^On- This is the same as saying that the naturally induced measure on 
T is normalized so that 9f{T) = 1. This normalization leads to the uniform 
distribution equation ([2]) in the form given in Thm.nbelowQ 

For W CH, 

A{W) :={u£L:u*£W}. 

A set il. C H is called a window if S° C O C S for some compact set 
T, C H which satisfies S° = S. We shall often deal with families of windows 



In particular situations, for instance if the internal space is a real space, one may 
prefer to normalize on the basis of what seems natural for 9h and thereby introduce a 
multiplicative factor into ((Sjl . This amounts to corresponding variations by multiplicative 
factors in the correlations and frequencies in which we are interested, but has no intrinsic 
importance to what we are discussing here. 
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based on a single S. A model set or cut and project set is a set of the form 

A{x, y) = A{x, y, := x + A{-y + , 

where i7 is a window and (x, y) G M'^ x H. The model set is called regular if 
the boundary of 17 (or equivalently of S) has measure with respect to the 
Haar measure of H. 

It is easy to see that for a fixed Q, A{x, y, Q) depends only on := (x, y) 
mod jC, so we may write A{^, Q) instead. 

The cut and project scheme S together with a compact window E is called 
irredundant if the equation i + S = S implies that t = 0. When dealing 
with model sets it is always possible to adjust the cut and project scheme 
(by factoring out a subgroup of H) to get an irredundant cut and project 
scheme which defines exactly the same family yl(^, O), S° C f2 C S, though 
now E is replaced by its image in the quotient of H [25l [H [16] . 

If S and S' are irredundant cut and project schemes for the same model 
set (or ones differing on sets of density 0) then their internal spaces H, H' 
are isomorphic topological groups by an isomorphism that induces an iso- 
morphism of the corresponding lattices in the obvious way. Thus one may 
speak of the irredundant cut and project scheme of a model set. The proof 
of this essentially follows from the construction of H given in ^ see \16\ [2] . 
We shall always assume that we are in the irredundant situation. 

Model sets A, regular or not, are Delone subsets of M'^, that is to say, there 
exist positive real numbers r,R so that the cubes Cr,CR of side lengths r 
and R, no matter where they are translated to in M'^, have at most one point 
of (respectively at least one point of) A. Since any model set A = A{x, y, J7) 
also satisfies A — A C A(0, — 0,) (sets of differences) and since ^2 — $7 is 
relatively compact, we see that A — A is also uniformly discrete, and by a 
similar arguments, all finite sets of sums and differences 

A lb • • • lb A (n terms) 

are also uniformly discrete. This is the Meyer property of model sets [19j . 
Model sets are uniformly distributed point sets: 

Theorem 1. [20j Let Vt d H he measurable and relatively compact. Then, 
assuming the normalization of measures assumed in ■ 

(2) hm -^card(yl(e, f7) n C^) = ^/f(17) 

for ^ G T, almost surely. If the boundary of has Haar measure then 
the result holds for all gT. □ 

Remarks: 

(i) When we write lim/j^oo £(^Cr) mean that we use a (any) discrete 
increasing sequence of positive real numbers {Rj} oo. 

(ii) In this paper we need to work with averaging sequences like the one 
in Thm. [TJ The results here apply for any van Hove sequence {An} 
of subsets of M'^ satisfying the condition that there is a constant 
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C > with £{An — An) < Ci{An). In the case of regular model sets, 
which are the main focus of this paper, all results quoted that depend 
on averaging (frequencies, correlations) are actually independent of 
which averaging sequence we use. This is a consequence of unique 
ergodicity. 

(iii) In particular, although all averaging results in this paper are written 
with averages over Cr, the centres need not be restricted to and 
cubes could be replaced by balls, etc. 



3. CORRELATIONS AND DIFFRACTION 

Let A = Q) be a regular model set. The frequency in A of a set of 
points {0, xi, . . . , Xn} C M.'^ is defined as 

freq({0,xi,...,a;„}) = lim card{?/ G Cr : y, y+xi, . . . , G A} . 

it— >oo f^y'^R) 

The frequency is the expected number of occurrences of the pattern per unit 
of vollumc in R'^. There is no need for the elements of {0, xi, . . . , x„} to be 
distinct, though repetitions can clearly be deleted. 

The following is a well known consequence of the uniform distribution of 
model sets: 

Proposition 1. Let Ti C H be non-empty and compact with S = S° and 
OnidT,) = 0. Let A = A{^,0,) be a regular model set with S° C ri C S. 
Then the frequency of {0,xi, ... , Xn} in A = A{(^, Q) exists and 

n n 

freq({0, xi, . . . , x„}) = 0H{n D f| (-x* + f^)) = (S n f| (-x* + S)) 

j=l 3=1 

for all ^ G T. In particular, the frequency o/ {0, xi, . . . , x„} in yl(^, fi) does 
not depend on ^ but only on the cut and project scheme S and the closure 
T, = ft of the window. 

Proof: First assume the simple case that A = A{Q,). Then y,y + 

xi, . . . ,y + Xn € A iff y*, y* + x^^, . . . , y* + X* G iff y* G W := ^ H 
f}j^i{—Xj + Q). The quantity we are looking for is then 

lim -^card(yl(I^) n Cr) = OniW) , 

which proves the first claim in this case. In the general setting, we are 
looking at x + A(—y + Q), but it is clear that these translations do not affect 
the outcome. Nor does replacing O by S, which only results in measure 
changes to the set VF. □ 
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The (n + l)-point correlation (n = 1,2,...) of a model set A (or more 
generally any locally finite subset of M*^) is the measure on (M*^)" defined by 

^ yi,...,yn,x£CRnA 



R^oo £(CbI 

yi,...yneA 

for all / € Cc((M'^)"). Here is simultaneous translation of all the variables 
by X. The simpler second sum is a result of the van Hove property of the 
averaging sequence {Cr}. Because model sets are Meyer sets, the sets of 
elements yj — x which make up the values of the arguments of / occuring in 
the sums, lie in the uniformly discrete set A — A. Hence for model sets we 
find that 

(3) 7a"^^^ = Y freq({0,j;i,...,x„})5(^^^,„^^^). 

a::i,...,x„eA— A 

In view of Prop. [H all the correlations of model sets exist and they depend 
only on the closure of the window. In other words, for a given compact 
window S which is the closure of its interior, all model sets A{^, with 
S° C C S have the same correlations of all orders. Thus we can collect 
the model sets into families J-{S, S) which all have the same correlations. 

The diffraction of a point set in is, by definition, the Fourier transform 

(2) 

of its 2-point correlation. In the case of model sets, this can be described 
explicitly. For the given cut and project scheme 5, see ([1]), there is a Fourier 
dual 5 of it: 



(4) U 
L° ^ f 

formed by taking the Fourier duals of the groups involved in S [19]. All 
of these groups, being duals of locally compact Abelian groups, are also 
locally compact and Abelian. Although is canonically isomorphic with 
W^, sometimes, for the sake of clarity, it is convenient to make the notational 
distinction, as we do here. 

Here tti and 1^2 again are projections. The dual of the compact group T 
is discrete, and 7ri||, is injective. The image of T is denoted by L° . There is 
again a mapping * : L° — > H. 

The diffraction of any regular model set defined by a window $7 is pure 

point p5| and is supported on a subset of L° of M'^; it is explicitly given by 

(5) 7f ({fe}) = |l5^*)l'. 
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The set of /c G M'^ at which the diffraction is not zero (and hence is 
actually positive) is the set of locations of the Bragg peaks of A (the Bragg 
peaks being the combined information of the location and intensity of the 
atomic part of the diffraction). 

We wish to prove that the 2- and 3-point correlations completely classify 
all such families J-{S, S); that is to say, given that a point set is a model set 
in W^, then the 2- and 3-point correlations determine the cut and project 
scheme S and the compact window T, up to translation. We begin with an 
abstract result on extending partially defined group characters. 

4. Extending partial characters 

Let G be a locally compact Abelian group and let U{1) denote the com- 
pact group which is the unit circle in C. The dual to G is the group G of 
all continuous characters, i.e. continuous homomorphisms of G into U{1). 

Suppose that E is a closed subset of G with no interior, and ^ E. Let 
D := G\E and L>(2) .- {(^^^ ^2) : /ci, fca, fci + fca S D}. 

Lemma 1. D^"^^ is dense in G x G. 

Proof: Suppose D^"^^ is not dense in G x G, i.e., there is a non-empty 
open set [/ X F C (G X G) \ L>(2), where U,V C G are open. Since is a 
closed subset of G with no interior, D nU, D CiV are nonempty open sets. 
For allu e DnU,v G DnV, we have {u,v) G {GxG)\D^'^\ i.e., u + v e E. 
Thus, D nU G —V + E. This is impossible since -D n C/ is an open set and 
E has no interior. □ 

Proposition 2. Let 

ip:D — > U{1) 
be a continuous mapping satisfying 

(p{s + t) = (p{s)(p{t) 

whenever s,t,s + t £ D. Then there is a unique character x ^ G with 
X\D = if. 

Proof: Let li be the uniformity on G defined by its structure as a topo- 
logical group: the basic entourages are the sets 

\J{V) ■.= {{x,y):GxG,x-yeV] 

where V runs through open neighbourhoods of G G. This uniformity 
also induces the relative uniformity on D by intersecting its entourages with 
D X D, and this relative uniformity induces the relative topology on D [15j, 
Ch.6. 

We claim that ip : D — > U{1) is uniformly continuous. We show that 
given any e > there is an entourage \J{V{e)) D {D x D) for which {s,t) G 
U(y(e)) r\{Dx D) implies that \ip{s) - ip{t)\ < e. 

In fact V{e) := {s £ D : |(/j(s) — 1| < e} works. This is an open subset of D 
containing and furthermore, {s,t) G \J{V{e))r\{Dx D) implies s — t£ l/(e) 
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and then s — t £ D and \(p{s — t) — 1\ < e. Using the basic relation satisfied 
by if, \(p{s) — (p{t)\ = \ip{s — t)ip{t) — ip{t)\ = \ip{s — t) — 1| < e, which what 
we wished to show. 

Since G is locally compact, it is complete (see Corollary 1 in Chapter 3.3, 
[5]). Since E has no interior, G is the closure of D. Since ip is uniformly 
continuous on D it extends uniquely to a uniformly continuous function 
X ■ G — > U{1). Then the mapping G x G — > U{1) defined by {x,y) ^ 
xix + y)x{x)^^xiy)~^ is continuous and is equal to 1 on all of the set 
By Lemma [J D^-^) is dense in G x G and so, by the continuity, this mapping 
is identically equal to 1. Thus x is a character. □ 

5. From correlations to model sets 

Model sets are an important modeling tool in the subject of quasicrystals. 
However, recognizing them is awkward because the cut and project schemes 
that underlie them are not obvious from the model sets themselves. 

Here we quickly go over a construction given in [3], that allows one to 
recreate an irredundant cut and project scheme for a model set A, using the 
2-point correlation measure 7^^^. 

We know that 7(2) is supported within the uniformly discrete set A — A, 
which is a Meyer set. Let L be the subgroup of M°' generated by the set 
A — A. For each e > define 

Pe:={te L: dens((t + A) A A) < e} , 

where A is the symmetric difference operator. We can use the P^, < e < 
2dens(A), as a neighbourhood basis of of a topology, called the autocor- 
relation topology, on L that makes it into a topological group. is called 
the set of the e— almost periods of A. 

It is a notable fact about model sets that the sets are relatively dense, 
a consequence of the uniform distribution, see Thm. 3 [3]. 

Now we define H to be the (Hausdorff) completion of L in the autocor- 
relation topology. Then a uniformly continuous homomorphism (p : L ^ H 
exists with the following properties: 

i) the image 'p{L) is dense in H; 

ii) the mapping 99 is an open mapping from L onto f{L), the latter 
with the induced topology of the completion H; 

iii) ker((/3) = closure of {0} in L. 

Moreover, since P^ is precompact in the AC topology, H is a locally compact 
Abelian group. 

Proposition 3. H is compactly generated. 

Proof: L is finitely generated ([.19j, Prop. 7.4), say with finite generating 
set J = —J. Then for any relatively compact open neighbourhood U of 
{0} in H, V := \Jj^j{U + (p{j)) is also open and relatively compact. Let 
Vn ■= V + ■ ■ ■ + V (n summands). Then IJ^i = H. In fact, if x £ H 
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then there exists y G {x — U) Ci (p{L) since (p{L) is dense in H. Writing y = 

\-jn) for some ji, ... ,jn ^ J,we have x = u+ip{ji)-\ ^fijn) G K 

for some u U. Thus V is compact and generates H. □ 
Now we define L := {(t,(p{t)) : t S L}. Then L is uniformly discrete 
and relatively dense in R*^ x H. Hence, L is a lattice of M'^ x H and S := 
{W^,H,L), along with the natural projection maps, is a cut and project 
scheme. We introduce the mapping * : L — > H as above. It is nothing 
other than (p. 

Theorem 2. Ti := A* C H is compact and is the closure of its interior. 
Furthermore, A = A{^}) for some 17 lying between T, and its interior. 

This result can be found in [2j Prop. 4 and Prop. 6. 

Although we see now that A is a model set from the cut and project 
scheme S, and we have constructed S from 7^^^, nonetheless, as we pointed 
out in the Introduction (see [IJ), j^'^^ does not contain enough information 
to determine the window. 

6. The role of the 3-point correlation 

Let A = A{Q) where is a window in H with boundary of measure 0. Let 
E = E{n) :={ke H : l^{k) = 0}, and D :=H\E. Since In is compactly 
supported and measurable, it has a continuous Fourier transform 1q. The 
set E is closed, so D is an open set. Moreover, since ln(0) = > 0, 
we have £ D. Measure alterations to 0, do not affect E and D. The 
relevance of D and E is immediate from ([5|). 

We are going to show that the 3-point correlation is decisive in charac- 
terizing a model set (amongst all other model sets) as long its window 0, in 
its irredundant cut and project scheme satisfies E{Q)° = 0. 

For real internal groups this is assured: 

Proposition 4. If H = then the sets E{Vt) have no interior points. 

This is a consequence of the Paley- Wiener theorem [24], Thm. IX. 12 since 
has compact closure. However, it is also easy to see directly: 
Proof: Let g be the real part of the function 1^- For k £ M"^, 

g{k) = I cos(27rA; • x)dx . 
Jn 

Writing out the Taylor expansion of cos around the origin and integrating 
term by term, it is easy to see that g has a Taylor expansion valid on all of 
M*^ and hence is an analytic function. It follows that its zero set in M*^ has 
no interior. Since In (A:) = only if g{k) = 0, we conclude that the set E 
has no interior. □ 

Let us return to the general situation with an internal group H. Let 
Q G H he a window in H and assume that E{Q) has no interior points. For 
notational simplicity let / := Iq. 
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Since D is open and G there is a ro > such that the cube C^o C D. 
We define functions I"" : — ^ M, n = 1, 2, . . . , by 

n 

i^''\wi,...,wn) ■.= e{nnf]{wj + n)), 

or equivalently, 

„ n 

(6) Z^^\wu...,wn)= \{hwj-t)f{t)dt, 

where f{w) := f{—w). 

Lemma 2. For n G N, T^") is uniquely determined by T^""*"^^- 

Proof: It is clear that X^") is a continuous function supported inside 
{Q, — Q)". Since A* is dense in it follows that (A — A)* is dense in — fi. 
Moreover, for x\, . . . ,Xn e A- A, 

T(-)(xt,...,rr*)=7i"+'^(-ari,...,-x„), 

which implies that T^"^ is uniquely determined by 7^*^^^- ^ 
Up to density changes, the compact set O is equally well described by 
its characteristic function / = 1q or by its Fourier transform /. Using the 

(2) (3) 

latter we shall show that is determined up to translation by ^\ and . 
Define ^0 on D by 

f{k) 



\m\ 

Then ^o(^) is a continuous function on D and |0o(fc)| = 1- Since /(O) 
£{n) > 0, ^o(O) = 1- A. simple computation shows that 



n 



(7) XH(fci, . . . , = n f{k,)fiY^ k, 

j=l j=l 

When n = 1 

T^^\k) ■= I f{k-t)f{t)dt. 

Thus J(^) is the convolution product of the function / and / and 

fW(fc) = />)/>) = |/>)p. 

When n = 2, 

(8) ^KkiM)=l{ki)l{k2)f{ki + k2). 

Let := {(A;i,/c2) : ^1,^2,^1 + A;2 G D}. By the definition of ro, 

Cro/2 X C'ro/2 C D^'^h On D^^) (Jggng 



\f{ki)mk2)\\m+k2)\ 
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By ©, for {ki,k2)eD^^\ 

(9) Mki + k2) = Mh)Mk2)i^^^Hki,k2). 

This imphes that the function cpo is a particular solution of the equation 

(10) Hh + k2) = Hhmk2)i^^^Hh,k2), 

where (/> is defined on D and {ki,k2) G D^'^\ We point out here that this 
equation is entirely determined by the functions X^^^T^^^ since the function 
■0(2) ig given by them. 

Equation (jlOh is related to the 'homogeneous' equation 

(11) ip{ki + k2)=ip{ki)ip{k2), 97(0) = 1, 

where f is defined on D and {ki, ^2) £ D^'^\ 

Let i;^ be an arbitrary solution of equation (llOh . Then i?i>/(/)o is a solution 
of equation (jlip and by Prop. [2] (with G there being replaced with H) it is 
the restriction to D of a character on H; that is, it is the restriction to H 

of an element of H ~ H. Hence it is determined by some element a & H 
through Xa ■ k ^ {k,a) £ U{1) and each solution of (|1U|) has the form 

(12) m-=Mk)Xaik). 
Finally, we get the main result of this section. 

Theorem 3. Let A = A(Q) be a regular model set and assume that E{^}) 
has no interior points (which is guaranteed if H is a real group). Then any 
model set with the same 2-point and 3-point correlation measures as A is, 
up to density 0, a translate of A. 

Note: All elements of the hull differ from translates of A by sets of density 
0, see [2j, Prop. 7. 

Proof: Let A' be a model set for some irredundant cut and project 
scheme S' and suppose that it has the same 2- and 3-point correlations 
as the model set A = A{il.) from the irredundant cut and project scheme 
S. The equality of the 2-point correlations shows that we can assume that 
S' = S. 

Recall that equation (jlOp is determined by X(^) and X(2). Since A' has the 
same 2- and 3-point correlations as A, the function (p for A' corresponding 
to the function 0o for A is another solution of (llOp . As we have already 
shown, this implies that (|12p holds for some a £ H. 

Let /' be the characteristic function of the window of A'. Then (p = f'/\f'\ 
and \f'\ = {lW)h = |/|. Thus 

f'{k) := f{k)xa{k). 
Taking the inverse Fourier transform on both sides of this equation we have 
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i.e. they are equal except possibly on a set of measure 0, which shows that 
up to density 0, A' is a translate of A. □ 

7. Counterexamples 

In this section we give examples of cut and project schemes and model 
sets from them for which the 2- and 3-correlations do not determine the 
window of the model set up to translation or inversion. For these sets the 
corresponding sets E{i}) have interior points. 

7.1. A periodic example. The most obvious way to make the sets E{^}) 
have interior points is to use a compact internal space, for then H is a 
discrete group. In the case that H is a finite group we are dealing with 
periodic structures, and for these it is long known that 2- and 3-correlations 
may fail to characterize the structure. The paper of Griinbaum and Moore 
[12j has useful introduction to the homometry problem for crystals (periodic 
structures) and some good ways of making examples. One such example, 
put here into the language of model sets, is the following: 

With N = 32, there are sets A,Bc {0, . . . , 31} which have the same 2- 
and 3-point correlations (when the point patterns are treated modulo N) 
[r2j §5.3, but which are not equivalent by 'rigid motions'. That is, they 
have the same pattern frequencies in Z/NZ for all 2- and 3-point patterns. 
Explicitly such a pair is given by the exponents of the expansions of the 
polynomials 

\ — 

pJx) ■= (1 - + X^)(l - X + - + x^) = 1 + x'^ + + + x^^ 

1 — X 

+^15 + ^17 + ^18 ^ ^19 ^ ^20 ^ ^21 ^ ^22 ^ ^26 + ^27 + ^29 + ^30 

1 — 

Pb{x) := (1 - x^ + x^){l - x"^ + x^ - x^ + x^) = 1 + X + x^ + x^ + 

1 — X 

+3;12 + + ^15 + ^18 ^ ^19 ^ ^20 ^ ^21 ^ ^22 ^ ^23 ^ ^27 ^ ^30 _ 

That is, A = {0, 7, 8, 9, 12, ... , 30} and B = {0, 1, 8, 9, 10, ... , 30}. Com- 
paring the gaps between successive positions in A and B along the line, 
namely (6, 0, 0, 2, . . . , 1, 0; 1) and (0, 6, 0, 0, . . . .3, 2; 1) where we wrap around 
modulo 32 at the semicolon, one sees directly that the gap patterns are nei- 
ther translationally equivalent or equivalent by reversal of direction. 

Now form the cut and project scheme 

M ^ M X Z/NZ ^ Z/iVZ 
U 

L ^ Z 

X < > (x,XAr) H-> Xat, 

X mod N . 



(13) 

where xv := 
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Figure 1. A plot over a complete cycle of the diffraction 
of the periodic model set determined by the set A of §7.1. 
The tick mark labels k stand for k/32 and the domain is 
shown for = 0, . . . , 32. Note the vanishing of the diffraction 
(extinctions) at all points 2k/32, k^O mod 16. 



Then for A and B, taken modulo N, we have A{A),A{B). These are 
periodic subsets of Z which arc located at the points of ^ + NZ and B + 
7VZ. To keep things straight below, we will denote them by Ap{A),Ap{B), 
indicating their periodic nature. The pattern frequencies for 2- and 3-point 
patterns of Ap{A) and Ap(B) are not altered when reduced mod A^, and 
hence are equal. Assuming normalization so that the density of Z is 1, the 
appropriate volume function on Z/NZ is yo\x/nz{S) = cajd{S)/N for all 
S C Z/NZ. 

This produces two (periodic) model sets which are not related by trans- 
lation or inversion of windows but which have the same 2- and 3-point 
correlations. The diffraction from either of the two model sets is the same. 
It is periodic of period 1 and supported on the set ^Z with values explicitly 
given by 

2 



-^PA(exp(27rzfc/32) 



at A;/32, A; G Z. (The label A can be replaced by B here). 

7.2. An aperiodic example. We can use this example to create distinct 

aperiodic model sets on the line which have the same 2- and 3-point corre- 
lations. The example below, based on the standard Fibonacci model sets, 
shows how this can be done. 

Let r := (l-|-\/5)/2, let L := Z+Zt, and let ' : L — > L be the conjugation 
mapping defined by r' := (1 — \/5)/2 . Form the standard Fibonacci cut 
and project scheme 

R 4^ M X R ^ M 

^^^^ L ^ C — > L 

X < — ^ {x, x') I— > x' . 
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Figure 2. A fragment the Fibonacci model set {. . . , — 1 — 
r, —1, 0, T, 1 + r, . . . } based on the window [—1, r^^) and the 
correspondingly thinned model sets using the congruence sets 
A and B of §7.1 above and below it respectively. These two 
thinned sets have the same second and third moments, al- 
though they are not translationally equivalent. The location 
of the origin is indicated by the black vertical line. 

Each X = u + VT £ L, u,v £ Z, is mapped to (x, x') = u{l, 1) + v{t, t') G C. 
We define a : L — > Z/NZ by a{x) = un, i.e. reduction of u modulo N. 
Let be a window in R and yl(H^) the corresponding model set. 
Now consider the combined cut and project scheme: 

M X (M X Z/NZ) ^ Rx Z/NZ 
U 

Cc — > Lx Z/NZ 

(x, (x', X* := (x', a(x)) . 

Notice here that Cc ■= {(x, (x', a(x))) : x G L}. With the subsets A, B 
above, we obtain the model sets Ac(W x A) and Ac{W x B), where the 
subscripts c stand for the combined cut and project scheme. In effect these 
consist of the points x of yl(VF) for which a{x) G A (resp. B). 

Now consider any pattern {0, x,y} C L for Ac(W x A). Write x* = (x', r), 
y* = {y',s) with r,s £ A. Up to the appropriate normalizing factor (see 
Thm. [T|), the frequency of the pattern in Ac{W x A), which is also the 
3-point correlation 'j^^\{x,y)), is 

freq({0,x,y}) = yoIc{{W x A) n {-x* + {W x A)) n {-y* + {W x A)) 
= volc{{W n (-x' + W)n i-y' + W)x 

An {-r + A)) n {-s + A)) 
= vo\K{iWni-x' + w)ni-y' + W)). 

Yoh/Nz{A n i-r + A)) n i-s + A)) . 

Replacing Ahy B gives the three point correlation at (x,y) for Ac(W x B). 
However 

woh/Nz{A n (-r + A) n {-s + A)) = voh/Nz{B n i-r + B) D {-s + B)) 

and so we have shown that Ac(W x A) and Ac(W x B) have the same 3- 
point correlations (and hence also equal 2-point correlations). Due to the 



(15) 
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different gap patterns produced by A and B the windows of the two model 
sets cannot be matched by translation or inversion. 

7.3. Aperiodic x periodic. Even simpler is to construct model sets of the 
form aperiodic x periodic with identical 2- and 3-point correlations Suppose 
that E := yl(VF) is a model set arising from the cut and project scheme ([I]). 
We take the direct product of this cut and project scheme and the one given 
at (lisp , and take as the windows the sets W x A and W x B. This leads to 
model sets S x A{A) and S x A{B). 

Each instance {z, z + zi, z + Z2} in S x A(A) of a pattern {0, zi, Z2} in 
M'^ X Z/A^Z is uniquely expressible as instances 

{x, X + xi, X + X2} X {a, a + ai,a + 02} 

from S and ^1(^4) of the product 

{0,xi,X2} X {0,01,02} 

of patterns from M*^ and Z/A^Z; and vice- versa. The frequency of {0, zi, ^2} 
in X TLjNTL is the product of the frequencies of {0, xi,X2} in S and 
{0,01,02} in A[A). 

The same goes when we use A{B\ and so S x yl(^) and S x A{B^ have 
the same 2- and 3-point correlations. 

7.4. Comments on real spaces x finite groups. Let F be a finite group 
(with the discrete topology) and F its dual (also finite). We are interested 
in the Fourier transforms of characteristic functions on subsets of spaces of 
the form M™' x and in conditions under which they may have zero sets 
with non-empty interiors. 

Give X F the product topology. Let C x F be a non-empty 
relatively compact set satisfying = 17° and write 

0= |J(Oa,o). 

aG-F 

Let A be the set of elements a £ F for which 0^ has a non-empty interior. 
Then 1^ : IT" x F — > C: 

ln{k,b) = ^ / e-^-'^-^^dx = ^^iZik) 

where Xa is the character on F defined by o G F. 

For each fixed b £ F, this is an analytic function of k on M™ ~ R™. Let 
Eh denote the vanishing set of 1q on M"^ x {b}. In order that In vanish 
on a set with interior, we require that some of the Eb have interior points. 
Being analytic on M"* x {6}, this requires 



Y,Xaib)lZik)^0 
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and hence that 

(16) ^^lf,„(A;)^0 

as a function on M'^. 

There are straightforward ways to make this happen. For instance, if all 
Qa, a G A, are equal then the requirement is simply that J2aeA^a-i^) ~ ^ 
for some b £ F, which is certainly possible. 

This is what is going on in the examples above. The equation X^^g^ Xa{b) = 
becomes 

where 6g{0,1,...31}. This equation holds whenever b is even and different 
from 0, since then the factor • oipA{x) evaluates to 0. Similarly 

for B. 

On the other hand, for the Penrose point sets arising from the Penrose 
rhombic tiling, the vanishing set E is indeed without any interior points. In 
this case the internal space is C x Z/5Z and the windows are the pentagons 
P, —tP, tP, —P where P is the convex hull of the 5th roots of 1 and the 
four listed windows are for congruence classes 1,2,3,4 mod 5 respectively 
\W]. Here it is easy to see that the zero condition ([16]) cannot possibly 
be satisfied. So rhombic Penrose point sets are determined by their 2- and 
3-point correlations. 

One final comment. One can raise the bar and ask about examples with 
equal 4-point or higher correlations. There are no examples of the type 
that we have produced here, because there are no finite cyclotomic sets like 
A, B above which are translationally inequivalent but have equal 4-point 
correlations, [12]. However, extending the investigation to weighted point 
sets, there are examples with equal second through fifth point correlations, 

8. Connections with point processes 

Start with the cut and project scheme ([1]). Let Vt <Z H he a. window (note 
the assumptions on windows, given in ^ for a model set A = A(ri) in M'^. 
Let S := Vt. Then A is a Delone set and there is a positive number r so that 
its distinct points are separated by at least the distance r. Let be the 
family of all discrete subsets of M*^ whose points are separated by at least 
the distance r. There is a uniformity on whose entourages are generated 
by the sets 

Ur^, := {($, e : ^ n Cb, C + and n C + C,] 

as R, e vary over M>o. 

The topology defined by this uniformity is a Fell topology, and so T>r is a 
compact Hausdorff space (^, Thm.l). The translation action of M'^ on 2?^ 
is continuous and the pair (Pr,M'^) is a topological dynamical system p3] . 
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The hull X (A) of A is the closure in this topology of the set of all translates 
of A. It is clearly a compact subset of "Dr and is also translation invariant. 
Thus {X{A),R<^) is a dynamical system in its own right. 

It is known [25] that X{A) is uniquely ergodic. Let /i denote the unique 
ergodic probability measure on it. There is a continuous M'^-equivariant 
mapping /3 : X{A) — > T, called the torus parametrization, from the hull 
into the group T of the cut and project scheme The inverse image 

in X(A) of any point ^ = {x, y) mod C is made up of sets of the form 
X + A{-y + 17') for some set Q.' satisfying 11° C 17' C S [25l [2] • If A is 
a regular model set, then so are all the elements of ^(A). Furthermore, 
assuming regularity, /i-almost surely the inverse image of ^ = (x, y) mod C 
consists of just one point and it is x + A{—y + S°) = x + A[—y + S). These 
elements are called non-singular elements of -'^(A). See [2], §5.3 for more 
details. 

Thus the point sets that make up X{A) are model sets and almost all of 
them are nan- singular. All of the point sets in X[A) have point correlations 
of all orders k, and these /c-point correlations are identical for all the elements 
of -'^(A). Thus we may speak of the fc-point correlation of ^(A). 

There is an obvious mapping N from the bounded measurable subsets of 
into L2(A:(A),/i) defined by 

N[A){^) = card(^>ny4). 

The mapping N can be construed as an ergodic uniformly discrete point 
process, with Jx{h) ^{^)dfJ' being the expectation for the number of points 
in A for a randomly chosen (according to the law fj,) point set $ G X{A). The 
function extends naturally to a M^-equivariant mapping N : Cc(M°') — >■ 
L2(X(A),/z)) defined by 

for all $ G X(A). 

Following [lOj, there has been increasing interest in studying the math- 
ematics of quasicyrstals by using ideas from the theory of point processes. 
Although quasicrystals (and model sets!) are, by their very nature, assumed 
to be highly ordered (correlated), and hence quite atypical from the point 
of view of the theory of point processes, nonetheless the techniques of this 
theory are applicable and quite effective. For a more detailed and compre- 
hensive exposition of this point of view see [7] . 

There is a one-one correspondence between the correlation measures of 
a typical point set and the moment measures of the corresponding point 
processes [6], §12.2. The correspondence between the 2-point correlation 
and the first reduced moment measure was utilized by Gouere [10] in his 
analysis of pure point diffraction and almost periodicity. In [7], we prove 
that an ergodic point process of uniformly discrete point sets is uniquely 
determined by its moments, or equivalently, all of its correlation measures. 
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This implies that solving the inverse problem for quasicrystals (i.e. deter- 
mining the law of the corresponding point process) is equivalent to knowing 
all of its A;-point correlations. In this language we have from Thm. [31 

Theorem 4. Let X{K) and X{h!) he point processes from regular model 
sets and suppose that their second and third moments are the same. Then 
their cut and project schemes may he identified. If T is the corresponding 
compact Abelian group and (3, (3' are corresponding torus mappings, then for 
each ^ G T the elements of (3~^{^) and I3'~^{^) all differ from one another 
at most on sets of density 0. 

Another approach to the way in which correlations link to the structure of 
the dynamics in the pure point case is given in [17J. Here the setting is a pure 
point ergodic uniformly discrete point process {X,fj,,N). The diffraction is 
almost always the same for the A £ X and it the Fourier transform of the 
first moment of the Palm measure of the process. The group M generated 
by the set S = —S oi positions of the Bragg peaks is the Fourier module of 
the process. There it is proved that if S + • • • + = M, where there are n 
summands, then the point process is determined by the 2-, 3-, . . . , {2n + 1)- 
point correlations. Thus if there are no extinctions {S = M) we need only 
the 2- and 3-point correlations. This result is not as strong as we have 
obtained here for real model sets with real internal spaces, since we have 
made no requirements on S. On the other hand, it goes well beyond model 
sets and works for general locally compact Abelian groups, and furthermore 
suggests a fundamental role for the extinctions in understanding how much 
of the dynamics is controlled by the diffraction. 
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